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We calculate the polarization function including only intra-Landau level correlation effects of 
striped Hall gas. Using the polarization function, the dielectric function, the dispersion of the 
plasmon and the correlation energy are computed in a random phase approximation (RPA) and 
generalized random phase approximation (GRPA). The plasmon becomes anisotropic and gapless 
owing to the anisotropy of the striped Hall gas and two dimensionality of the quantum Hall system. 
The plasmon approximately agrees with the phonon derived before by the single mode approxima- 
tion. The (G)RPA correlation energy is compared with other numerical calculations. 

PACS numbers: 73.43.Lp 

I. INTRODUCTION 

At the half- filled third and higher Landau level (LL), the anisotropic longitudinal resistance has been observed 
in ultrahigh mobility samples^^ Collective excitations and related physical quantities of the anisotropic state give 
important information to clarify this state, but have not been studied well experimentally and theoretically. In this 
paper, we study collective excitations by calculating a polarization function in intra-LL. 

The Hartree-Fock approximation (HFA) at a half-filled Zth LL in I — > oo limit predicts that a unidirectional charge 
density wave (UCDW) forms in the two-dimensional (2D) electron systemA^ Within the intra-LL, the HFA at the 
half-filled ith LL (I > 2) shows that an anisotropic Fermi surface forms in the UCDW which explains the anisotropic 
longitudinal resistance and exotic responseA^ Collective modes for the UCDW are studied based on an edge current 
picture^ and a single mode approximation (SMA) J£ We call the UCDW the striped Hall gas. Another solution in the 
HFA is a modulated striped state which is highly anisotropic charge density wave (ACDW) having an energy gap and 
no Fermi surface. Collective modes for ACDW is derived by the time dependent HFA (TDHFA)£ In this paper, we 
study the striped Hall gas by using a random phase approximation (RPA), in which bubble diagrams are summed up, 
and a generalized RPA (GRPA), in which bubble and ladder diagrams are summed up. The spectrum of collective 
modes and correlation energy of the striped Hall gas are compared with the results of SMA and ACDW. 

The electron gas system in the absence of a magnetic field has been studied successfully by using diagram techniques 
in many-particle physics. Diagram technique enables us to compute corrections to HFA systematically^ In a 2D 
electron system in a strong magnetic field, on the other hand, the electron kinetic energy is frozen to one LL and it 
seems difficult to analyze the system by using diagram techniques. In the von Neumann lattice (vNL) formalism, 12 
however, it is possible to apply a systematic diagram technique to the 2D electron system under a magnetic field. In 
the vNL formalism, the 2D electron system under a homogeneous magnetic field is represented as a 2D lattice system 
and a momenta is defined in the magnetic Brillouin zone (MBZ). The electron kinetic energy is induced in HFA within 
intra-LL and depends on the momentum. Hence the free propagator is defined by using the kinetic energy and a 
diagram technique is applied systematically. 

In this paper, we investigate quantum fluctuation effects of the striped Hall gas below the cyclotron energy scale. 
From the polarization function including only intra-LL effects at one-loop order, we obtain the dielectric function, the 
plasmon and the correlation energy in the (G)RPA. Isotropic screening effects in the dielectric function at higher LL 
have been estimated^ The intra-LL screening effects in the dielectric function is highly anisotropic for the striped 
Hall gas. The plasma frequency which is obtained by the zero of the dielectric function in (G)RPA, is found to be 
an anisotropic gapless mode. In the RPA, the energy of plasmon is larger than the particle-hole excitation as a usual 
plasmon in the electron gas without a magnetic field. In the GRPA, on the other hand, the energy of the plasmon 
becomes smaller than the particle-hole excitation. The plasmon in the GRPA approaches to the phonon derived by 
the SMA as including a long-range component of the Coulomb interaction. The correlation energy in the (G)RPA 
substantially reduces the total energy. 

This paper is organized as follows. The striped Hall gas with the Fermi surface of a strip shape is introduced in 
Sec. [n] In Sec. IIIII we calculate the polarization function at the one-loop order and derive the dielectric function, 
the plasmon and the correlation energy in the (G)RPA. Summary is given in Sec. I1VI Appendix [X] gives the LL- 
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projected HF approximation in the vNL formalism. The explicit one-electron energy of the striped state is given in 
Appendix [BJ In Appendix lei the Feynman rule for the diagram techniques are presented. In Appendix [DJ a duality 
relation between the direct term and the exchange term is shown. 



II. STRIPED HALL GAS WITH AN ANISOTROPIC FERMI SURFACE 



In this section, we derive the striped Hall gas at the half-filled higher LL in the vNL formalism^A^ In the vNL 
formalism, the 2D electron system in a perpendicular uniform magnetic field is transformed into a 2D lattice system 
and the momentum becomes good quantum number of the one-electron state. In the HFA, the kinetic energy is 
induced by the interaction term and perturbative calculation is easily carried out by using the diagram technique. 

Let us consider a 2D electron system in a perpendicular uniform magnetic field B. The total Hamiltonian H of the 
system is written as H = Hq + H i , 



//,, / V f (r,f) iP+ n eA) \ (r,t)d 2 r, 
2 m 



Hi = - 
2 



p(r, t)V{v - r')p{r', t) : d 2 rd 2 r' 
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where p a = —ihd a , d x A y — d y A x — B, V — q 2 /\r\, q 2 — e 2 /Ane (e is a background dielectric constant). ip(r,t) is the 
electron field, and p(r,t) = i/j'(r,t)^j(r,t). The symbol : means the normal ordering for the vacuum of one electron 
field. In the following calculation, the units are set as h = c = 1. Hq is the free Hamiltonian, which is quenched in a 
LL. Hi is the Coulomb interaction. Since we want to consider only intra-LL effects, the free Hamiltonian is omitted. 
We ignore the spin degree of freedom. The electron field is expanded by the momentum state |// (g> /3 p ) in the vNL 
formalism 1 - as 
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where 6;(p, t) is the anti-commuting annihilation operator with momentum p defined in the MBZ. The base function 
depends on two momenta p x , p y symmetrically. The annihilation operator &;(p, t) obeys a twisted periodic boundary 
condition bi(p — 27rn,t) = e~™( nx+ny * >+lnyP!C bi(p, t), where n x , n y are integers. The momentum state is the Fourier 
transform of the Wannier basis of vNL which localized at r = a(r s m,n/r s ), where n, m are integers. Here a = 
^j2irh/eB, and r s is a vNL asymmetry parameter. In this paper, we set a to be 1. The MBZ means \p a \ < tt. The t 
dependence of the operators is omitted for simple notation. The Fourier transformed current operator j M (k) is written 
in the vNL formalism as 



J"(k) = 
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where k = (r s k x , k y /r s ) and £ is the relative coordinate of the electron and = [1, —(eB/m)r], (eB/m)£], p =0, 1, 
2. In order to project the operators into the Zth LL, we take only Zth LL index in Eq. J2J and write it as J* (k). We 
define gf(fc) = (fi\^{v^, e~ lk '^}|/;}. For the projected density operator /0*(k) = j°(k), we use the short-hand notation 

fi(k) = g®(k) and fi(k) = Li(k 2 /4n)e~ k / 87r where Li is the Laguerre polynomial. In the LL projected space, the 
Fourier transformed Coulomb interaction is modified as Vi(k) = [/;(fc)] 2 27rg 2 jk for the vNL operator 6j(p). We apply 
the HFA to the Coulomb interaction within the Zth LL, and get the one-electron spectrum satisfying a self-consistency 
explicitly (see Appendix . 

There are two conserved charges Qx and Qy corresponding to the magnetic translation in — y direction and x 
direction, respectively, 



Qx 
Qy 



. dp* (k) 



dk x 
.9p*(k) 
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k=0 



k=0 



which satisfy [H, Qx] — [H, Qy] = 0>i& Note that Qx and Qy correspond to the magnetic translation in p x direction 
and p y direction in the momentum space, respectively. We give a self-consistent mean field solution which is uniform 
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in the y direction and periodic in the x direction, that is the striped Hall gas. This state is given as 

|HF)=iV 1 J] &[(p)|0), (5) 

|Pa:|<1',b H |<1'/2 

where |0) is the vacuum for bi and N% is a normalization factor. This striped state satisfies the self-consistency 
Eq. (|A3j) at the half- filled higher LL£^ The corresponding one-electron energy has the anisotropic Fermi surface 
which is parallel to the p x axis. The explicit form of ep F (p y ) is given by Eq. l|B2jl in Appendix FBI The density of this 
state (HF|p(r)|HF) is uniform in y direction and periodic in x direction with a period r s m^ Using Eq. (j3J), we can 
show 



(HF|[Q x ,p,(k)]|HF) = 0, (6) 
<HF|[Q^,p*(k)]|HF) ^ Ofor (k x ,k v ) = (27rn/r.,0), 

where n is integer. Therefore the magnetic translational symmetry in x direction or p y direction is spontaneously 
broken. 

The one-electron energy has an anisotropic energy gap in p x direction. The Fermi velocity, then, is in the y 
direction in coordinate space. The orthogonality of the Fermi surface in the momentum space and the density profile 
in the coordinate space is reminiscent of the Hall effect. The Hartree-Fock (HF) energy per particle is calculated as 
i? ; HF (r s ) = (HFjiJi |HF)/iV. We determine the optimal value r s = r" lm which corresponds to the stripe period by 
minimizing £^ HF (r s ). At the half-filled I = 2 LL, the optimal value is r™ in = 2.474£ 



III. INTRA-LL POLARIZATION EFFECT 



In the Ith LL Hilbert space, the free Hamiltonian Hq is quenched and only the interaction Hamiltonian Hi remains. 
Since there exists no bare kinetic term, it seems difficult to deal with the Coulomb interaction as a perturbative term, 
naively. In the HFA, however, the kinetic term of electrons is induced by the effective direct and exchange interaction. 
Hence, the effective kinetic term appears in the HFA and the quantum fluctuation around the HF ground state is 
caused by the residual interaction term. The Coulomb interaction Hi is divided into two terms 

Hi = ff HF + (Hi ff HF ) = i/ HF + \ ° Mm ( k )M-V ° . (7) 

Here, £f HF is given by Eq. i|Al(l in Appendix lAl Operators between the symbol ° are normal ordered for the HF 
vacuum. We study the quantum fluctuation for the HF state using the polarization function in the (G)RPA. In the 
vNL formalism, we use the Feynman diagram technique in the vNL formalism which is presented in Appendix [C] 
The dielectric function and excitation spectra of the striped Hall gas are given by the polarization function in the 
(G)RPA. 



A. One-loop polarization function 



First let us study the polarization function in the one-loop order. The current-current correlation function, which 
is a response of the external electromagnetic field, in one LL is defined in the Hciscnbcrg picture as 
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t^ob^l^o) and TS is a total time times a total area of a 2D electron system. The one-loop 



current-current correlation function shown in Fig. ^ is calculated as 
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FIG. 1: Feynman diagram for the current- current correlation function. 



Gp°L is the free electron propagator given by Eq. (|C5I) . fip is the Fermi energy. We define the one- loop polarization 
function in the Ith LL by 7ri_i oop (k, uj) — 7r 00 (k, uj) / (g® (k)) 2 . On the basis that et is a even function, we can show 
that 7ri_i 00 p(-k, uj) = 7Ti_ loop (k, uj) and 7Ti_ loop (k, -uj) - 
part of the one-loop polarization function are as follows: 



7Ti-ioop(k, ui). The explicit form of the real and imaginary 



ReTT 1 _ loop (ky,u}) = 2 lim / — - 
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The real part is given by the principle integral. The delta function of the imaginary part means the energy conservation 
for the particle-hole excitation. The ^-independence is the result of the anisotropy of the striped Hall gas. The 
numerical results of 7ri_i oop (fc y , to) are given in Fig. [5] The uj region in Fig. [5] where Im.iri-i oop (ky,uj) takes finite 
value corresponds to the particlc-hole excitation region. At the lower boundary of the particle-hole excitation region, 
Im7ri_i 00 p(/cj,, uj) approaches zero. At the upper boundary of the particle-hole excitation region, Ini7ri_i oop (fcj,, uj) 
becomes infinite negatively. In contrast to the ordinary 2D electron gas system in which the Fermi surface is sphere 
shape, the particlc-holc excitation region has a gap for a finite k y due to one-dimensional nature of the striped Hall 
gas or the strip shape of the Fermi surface. 



B. Random phase approximation (bubble and ladder diagram) 

The quantum fluctuation beyond the one-loop order is calculated in the RPA which is the summation of the 
geometric series of the one-loop polarization function as shown in Fig. [3] Under a homogeneous magnetic field, we 
can also sum up bubble and ladder diagrams as shown in Fig.0J which is called the generalized RPA (GRPA). 15 As 
an example, the two- loop diagram is calculated in Appendix [D] 
The RPA polarization function is written as 

7TRPA(k, W ) = 4°J(k,uj)+Y / 4?A^Lo)W n (k)n^(k,uj) 

n 

+ 2^ ) (k,w)W'„(k)7rW (k jW )W TO (k)7r^(k,w) + • • • 

urn 

= 5^7r^ ) (k,a;)[l-W(A)7r( )(k,a;)]- 1 , (12) 

where we define 

Tr^fk u) = -i f°*E / d2p r (0) G (0) e ipx < n - m ' (13) 

W n (k) = V;(k + 27rn), and a x b = a x b y — a y b x . Here n = (n x /r s ,r s n y ). Since G p p is independent of p x for the 
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FIG. 2: ^-dependence of the polarization function 7ri_i oop (fcj / , w) at fc H = 7r/4. The unit of k y is r s /a. Im7ri_i oop vanishes at 
small uj region due to the strip shape of the Fermi surface. The unit of uo is q 2 /{ah). 




FIG. 3: Bubble diagram summation (RPA) 

striped Hall gas, we can integrate Eq. over p x , and then 

ir^Uk w ) = -j / ^ G (0) G (0) e -*P.(«--"»-)A n (U) 

iTnn is independent of n and equivalent to the one-loop polarization function 7ri_i 00 p(fc a , ui). Here, = Sn^oo 

appears as the result of dividing the infinite momentum integral region bytheMBZ. [l - W(k)irW (k, w)] ' _ means 

the (nO)-element of the inverse matrix of [l — W(k)ir^ (k, cj)1 . In the (G)RPA, the fluctuation for the transverse 
to the stripe are included through the argument of exp in Eq. {Sll - When the integer n x increases, the effects of the 
Coulomb interaction between different stripes increase because p y direction corresponds to the x direction. 

A peculiar property of a 2D system in a magnetic field is that the ladder diagram take a similar form with the 
bubble diagram in one LL^ In the presence of a magnetic field, there exists a duality relation between the direct 
term and the exchange term (see Appendix [D|. Owing to this property, bubble and ladder diagrams are able to be 
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FIG. 4: Bubble + Ladder diagram summation (GRPA). The dotted and dashed line means W and — W, respectively. The 
wavy line is W c s- 



summed up to the infinite order. The intra-LL polarization function in the GRPA is written as 
^ GRPA (k,c) = 4 ) (k,c,)+^ 7 r^(k, W ){^^fc)-^ l (fc)}7ri ) (k,c,) 

n 

+ J24n(.Ku){W n (k) - W n {k)}^± (k,uj){W m (k) - W m (k)}^ Q (k,uj) + ■■■ 

nm 

= ^4° n ) (k ) o;)[l-Ty eff (fc)7r( )(k )W )]; o 1 , (15) 

n 

where we define Wf{k) = W n (k) - W n (k). Here, W n (k) = Vi{^ + n y , ^ + n x ), and Vi(a) = J ^Vi(/3)e^ 3 . The 

negative sign in front of W n {k) is due to the exchange between the electron and hole in the ladder diagram. 

The effective interaction W§ s (k) is shown in Fig. at the I = 2 LL. The effective interaction W§ s (k) is positive 
for small k and negative for large k. Thus the bubble diagram is dominant in small k region and ladder diagram 
is dominant in large k region. As shown later, the negative contribution of the effective interaction in the ladder 
diagram changes the behavior of the dielectric function and plasmon drastically. 



C. Dielectric function 



We are able to obtain useful physical information of the system through the dielectric function. The RPA dielectric 
function is the denominator of the polarization function in the RPA. In the RPA, the pole of the polarization function, 
which is equivalent to the zero of the RPA dielectric constant, provides a spectrum of the plasmon excitation. 

In the RPA, the dielectric function is defined by 

e RPA (k,cu) = 1 - WofkjTri-ioopfo.w), (16) 

and in the GRPA, it is defined by 

e GRPA (k uj) = 1 - WofeW-ioopfe^) 

1 + Wo(k)7Tl_ lo0 p(fc y , uj) 

Here, we consider only n — term in Eq. (|15fl for the first approximation. The numerical results of the dielectric 
function at the RPA are shown in Fig. at typical k-value (k x ,k y ) = (ir/A, 1.5). The imaginary part e RPA of 
e RPA (k, uj) always takes a positive value, and the real part e RPA of e RPA (k, uj) has two zeros. e RPA (k, uj) is finite at 
the particle-hole excitation range. 

In the GRPA, the numerical results of e GRPA are shown in Fig.0 As shown in Fig. El the effective potential W§ s (k) 
has a negative region in contrast to Wo(fc), which causes the drastic change of the dielectric function behavior. Thus, 
the point where both e RPA and e RPA are zero, which corresponds to the pole of the plasmon, moves over the particle- 
hole range in GRPA (compare Fig. with Fig.0). 
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D. Plasmon 



The pole of the (G)RPA polarization function gives the excitation mode associated with the charge fluctuation, that 
is plasmon. The pole of the (G)RPA polarization function 7T( G ) RP a is zero of e( G ) RPA (k, uj). The plasmon appears at 

the outside range of the particle-hole pair regime where e^ RPA (k, uj) takes finite values. 

First, we see the case of considering only n = term of Eq. (|12[) and Eq. Q15[l. The plasma frequency given by 
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FIG. 7: ^-dependence of e GRPA which includes both bubble and ladder diagrams for {-^-g, 1 ' 5r * )■ The unit of uj is q 2 /{ah). 



solving the pole of 7T(g)rpa is shown in Fig. IHland Fig. El At k x ^ 0, the plasma frequency always approaches zero as 
k y — > 0. On the other hand, at fc a ^ the plasma frequency remains a finite value at k x = 0. The difference of the 
plasmon behavior between k x and k y direction is the result of the spontaneous breaking of the magnetic translation 
and rotation symmetry of the striped Hall gas. For the long wavelength limit and £ Py +k y — £p y *C u> p , the plasma 
frequency cu p rises like \k y \^/ ' — In fcyl/lkl 1 / 2 for taking only n = term. The origin of the square root behavior is 
the Coulomb interaction 1/r in two dimensional space. The logarithmic correction is caused by the divergent Fermi 
velocity due to the Coulomb interaction. 

For n — case, in the GRPA, the k y dependence of the plasma frequency separates two region. The one region is 
lower than the particle- hole pair excitation region, where W§ s (k) is negative. The another region is higher than the 
particle- hole pair excitation region, where W§ s (k) is positive. At the positive interaction region where bubble diagrams 
are dominant, the plasma frequency appears above the particle- hole excitation and collapses into the particle- hole 
excitation. On the other hand, at the negative interaction region where ladder diagrams are dominant, it appears 
below the particle-hole excitation. This negative interaction dominant state is considered as a low energy bound state 
due to the effective attractive interaction W§ s (k). At the long wave length range, the plasma frequency behaves as 
the same as the case of n = RPA. 

Next we include finite n terms of Eq. i|12fl and Eq. i|15|) . In the RPA, the plasma frequency is slightly larger than 
n — case at large k y region (see Fig. 0. On the other hand, in the GRPA, the plasma frequency is smaller than 
n = case for wide k y region (see Fig. |5J|, and approaches to the phonon frequency associated with the density 
fluctuation derived in the SMAiiS In the striped Hall gas, the phonon is a Nambu-Goldstone mode associated with 
the spontaneous symmetry breaking for the magnetic translation due to Qy in Eq. Increasing n in Eq. I|15|) 

means that the polarization function includes effects of nth MBZ increasingly. For the large n value, the argument 
of W° s (k) is in the large k range as seen in Fig. [SJ Since the charge density is the same as the electron density, 
it is reasonable that the plasmon associated with the charge fluctuation is the same as the phonon associated with 
the density fluctuation. For the small k y range, the convergence of the numerical calculation for large n is not good 
because of the singular behavior of the polarization function near k y = 0. 
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FIG. 8: fc^-dependence of plasma frequency uj p (k x , k y ) for the RPA in the n — case. The unit of uj p (k x , k y ) and k x is q 2 / (ah) 
and 1/a, respectively. 



E. Correlation energy 

We calculate the (G)RPA correlation energy in this section. As is well known, the correlation energy is given by 
the virtual coupling constant A integration^ 



^total _ £jHF 



1 I" 00 d 2 k ~ 

dX 2 J „ (2^ y(k)(i?(A)l ° Wfl-k) o \E(X)), 



(18) 



where p(k) = j°(k), V'(k) = is the Fourier transformed Coulomb interaction, and l-E'(A)) is the ground state for 
the system with the virtual coupling constant A. The second term of Eq. i|18fl is the correlation energy. By replacing 
p(k) with /5*(k) and V(k) with Vz(k) in Eq. H18[l. the LL projected correlation energy is represented by a vNL basis: 
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The integer n is caused by dividing the fc-integral region into summation of one MBZ. The RPA correlation energy 
is the sum of the chain diagram of bubble as shown in Fig. I1UI which is derived by the perturbative calculation about 
the virtual coupling constant A. The density correlation function part in Eq. (|19|) is calculated as 

d 2 Pi d 2 p 2 
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(20) 
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FIG. 9: fcy-dependence of the plasma frequency u v (k x — 7r/4, k y ) for the RPA and the GRPA with the particle-hole excitation 

region and the phonon frequency. The number means the dimension of the matrix |l — W aS {k)n^ (k, w)l . The dashed line 

is the plasmon in the RPA, and the dotted line is the one in the GRPA. The solid line is the phonon derived in the SMA^& The 
behavior of u> p approaches to the phonon one in the SMA when the matrix dimension increases. The particle-hole excitation 
region is unchanged and remains in the (G)RPA. The unit of u> p (k x , k y ), k x and k y is q 2 /(ah), l/(r s a) and r s /a, respectively. 
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FIG. 10: The chain diagram for RPA correlation energy. 
Here itnm is given in Eq. (|13|1 . We take only the diagonal matrix elements contribution and approximate Eq. (|20|1 as 

i E(^(*)Wn(*)) a - * - ■ (21) 

„ 1 - ir^ikjXWnik) 

In the summation about n, the contribution of n > 1 terms are negligible because of the Gaussian factor in Vi(k). 
So we consider only the n = term. By substituting Eq. (|2TJ> into Eq. l(T^|) . the integral about A gives the RPA 
correlation energy E^p" A as 

1 f°° iit,i f d 2 k 

EnPA = - 5 / 7T / 7^4 {log(6 RPA (k, w )) +n 1 - loop (k y ,u,)W (k)} . (22) 

2 J-oo 2vr J MBZ [2tt) j 

The numerical estimates of _E HF , the real part of E^p\ and the total energy _B total per particle at the half-filled I = 2 
LL are obtained as follows: 



E HF = 


-0.7706, 


TpCOTT 

^RPA — 


-0.0341, 


^total 


-0.8047, 



(23) 

where the energy unit is ^- . The total energy is lowered by the RPA correlation energy significantly. 
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FIG. 11: The chain diagram for GRPA correlation energy. 

GRPA correlation energy is the sum of the chain diagram of bubble and ladder as shown in Fig. El Corresponding 
to Eq. in RPA, Eq. is approximated in GRPA as 

« EOr&WwWWf (k)- 



l-7T^(fc)A^ ff (k)' 

Considering n = term, the GRPA correlation energy i?GRPA i s written by iri_i oop (k y , oj) as 

E C G °£>A = ~ fj^ ^^\(^^ (24) 

The numerical estimates of the real part of E(?£p A and the total energy £; total per particle at the half-filled / = 2 LL 
are obtained as follows: 

£grpa = -0.0385, (25) 
^total = _ .8090. 

GRPA correlation energy is slightly lower than RPA correlation energy. 

Yoshioka calculated the corresponding HF energy for the ACDW stated The ACDW state has anisotropic energy 
gaps in p x and p^-direction. Our striped Hall gas, on the other hand, has an energy gap only in p^-direction. The 
numerical value of Yoshioka's HF energy is E Y = —0.7763. Our total energy E total of the ground state is smaller 
than E Y , hence our striped Hall gas including the quantum fluctuation at the (G)RPA level is more stable than the 
ACDW. On the other hand, Shibata and Yoshioka studied the ground state phase of 2D electrons in I = 2 LL by a 
density matrix renormalization group (DMRG) method, which is a numerical calculation of a small system improved 
by the renormalization group method^ The DMRG results seem to predict the striped Hall gas at the half-filled 
higher LL. The total energy given by the DMRG is E D = —0.796 ± 0.004iS which is smaller than E Y , whereas it is 
close to £: total . 

The one-electron energy of the ACDW state has a gap in both p y and p x direction, whose value is about 1 Ki£ Hence, 
this state is insulator in x and y direction and the gap structure causes the quantization of the Hall conductance. 
However, experiments show the huge anisotropic resistivity and the Hall conductance is not quantized at several 
mKii*2 Since the striped Hall gas has the anisotropic Fermi surface, the x direction is insulator whose gap energy is 
the cyclotron energy, and the y direction is metal in which the electron gas state realizes. Therefore the striped Hall 
gas is more consistent with experiments than ACDW. Moreover the comparison of the correlation energy with the 
ACDW and results of DMRG seems to support our striped Hall gas. 
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IV. SUMMARY 



Using the one-loop polarization function which includes only intra-LL effects, the dielectric function, the plasma 
frequency, and correlation energy are calculated in the (G)RPA for the striped Hall gas. The characteristic feature of 
the plasma frequency is anisotropic gapless behavior. The anisotropy is due to the spontaneous breaking of rotational 
symmetry and the gapless feature comes from two-dimensionality of the system. The anisotropic plasma frequency 
will be observed by some experiments: e.g. a surface acoustic wave. The numerical result of GRPA plasmon suggests 
that the plasmon in the striped Hall gas is the same as phonon in the striped Hall gas which is the Nambu-Goldstone 
mode due to the the spontaneous breaking of translational symmetry. In contrast to the quantum Hall smectic 9 
derived by edge current picture and TDHFA applying to the ACDW, this excitation state reflects the striped Hall gas 
state. It is shown that the quantum fluctuation effect for the striped Hall gas substantially reduces the total energy 
in (G)RPA. This means that the quantum fluctuation plays a important role in the striped Hall gas. The quantum 
Hall gas properties strongly depends on the electron self-energy with the anisotropic Fermi surface. The treatment for 
quantum fluctuation effects to the electron self-energy is beyond the scope of the present paper, and is very interesting 
as a future problem. 
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APPENDIX A: HARTREE-FOCK APPROXIMATION 



In the intra-LL HFA, we reduce the Coulomb interaction term to the kinetic term by using a mean field 
U(p) = (HF|b|'(p)6i(p)|HF) where |HF) is a many-particle state satisfying a self-consistency equation. The inter- 
action Hamiltonian projected into the Ith LL, 

is approximated by the HF Hamiltonian 

% ?^r F (P - p'Mp'W(p)Mp) - ^ / r f^ r f^ (ph HF (p - p'Mp')- (Ai) 



A 7mbz(2 
Here, we define 



.HF 



(p -p')=J2 {Vi(27rn)e l (P-P') xn - V^irh + p' - p)} . (A2) 



v* is a filling factor for the highest LL, and N is a number of electrons at the highest LL. In this paper, v* is set to 
1/2. The states below I — 1 th LL are fully occupied. The first term in Eq. (|A2(1 is the Hartree potential and the 
second one is the Fock potential. The uniform positive background charge cancels the n = term in the Hartree 
term. Since the state at the Ith LL is occupied by electrons whose energy is below a Fermi energy /ip, the mean field 
is written as ?7(p) = ■^■9(hf — e j IF (p))< Hence, the self-consistency equation of one-electron energy e[ iF (p) reads 



™ = / £^«r (p - p') 9(,, F ef F (p>) ). (A3) 
The one-electron energy has a periodic structure e[ iF (p) = ef F (p + 27rn) owing to v ; HF (p) = «; HF (p + 2™). 
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APPENDIX B: ONE-ELECTRON SPECTRUM 

The one-electron spectrum is given by the next explicit relation ep F (pj,) = (p y ) + ef(p y ), where 



2j -Y, v ' 

oo 

E 



2 °° 



n=0 



2?r(2n+ 1) 



(-l) n cos[(2n+l)p. y ] 



(2n+l) 



* dfc. 



n—~oo 2 



y 

2tt 



.2tt 



\Jk 2 x + r 2 s {k y - Py - 27m) 2 



(Bl) 
(B2) 



APPENDIX C: FEYNMAN RULE 

In the following calculation, the interaction picture is applied in perturbation theory. In the Heisenberg picture, 
the time-dependence of an operator £> H (i) is defined as H (i) = e lHlt Oe~ lHlt . 
The electron Green' function is defined by 



(*o|T6£ )tl &£U|* > = iG(pi,ti - ia) ^(27r) 2 ^ 2 )( Pl - p 2 - 2™)^' 



(CI) 



where |^q) is the exact ground state of H\, and <j>(p, n) = n(n x + n y ) — n y p x . T means the time ordering. We write 
bi(p,t) as 6 Pj t and omit its LL index. The Heisenberg picture is changed to the interaction picture by using the 
relation bg t = S(0,t)b k . t S(t,0), where S(t,t') = e iH HF t e -iH lte iH™t' e -W x t' Using in Eq the Green 

function is written as a familiar form: 



G(k,t-t') 



1 > 

71=0 



00 (HF|T6 k X. t ^(^i)---^»)|HF) 
-oo *" (HF|5(+oo,-oo)|HF> : 



where the interaction V(t) is the residual interaction in the interaction picture 



v(t) = - 2 



' d 2 k 
o(2^ 



°£(k,t)iK*)p.(-Mo 



(C2) 



(C3) 



The lowest order Green function is obtained as 

G (0) (k, t - *') = -i U{t - t')9(e k - /i F )e~ l(t - t,)(ek -^ F) - 9(t' - t)9(fi F - e k )e" j( ^ t ' )(ek -' 1 f) } . 



(C4) 



Considering only intra LL effects, we omit the LL index and take the short notation of ep F (k) as e k . The Fourier 
transformation of the free propagator reads 



G 



(0) 



dt e^*G (0) (k,i) 



#(mf ^ £k) 



ui — e k + a*f + i$ uj — ek + [ip — i5 



(C5) 



In the vNL formalism, it is transparent way to represent the infinite k-integral of Eq. (|C3|I as the summation of 
one fundamental MBZ. So the residual interaction is written as 



V(t) = 



d 2 k d 2 pi d 2 p2 
MBZ (27r)2 (2ir) 2 (2tt)2 



Vl{Pl,K,P2)e o ° t pi+k °t,Pl t,p 2 °t,p 2+k oi 



where we define 



VS(pi,fc,j>2) = J2 V l (k + 27rh)e i ^-^ xa . 



(C6) 



(C7) 
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In stead of the infinite A:- integral region, the infinite summation appears. The local interaction Vj(fc) in momentum 
space is replaced by the non-local one Vi{p\, k,p 2 ) including the phase factor in the density operator. For perturbative 
calculations, we use the Wick's theorem and obtain an n-point correlation function for the interaction V(t). 

In the following we present the Feynman diagram rule in the momentum space for the perturbative residual inter- 
action Eq. (|C6|I . The phase factor due to the magnetic field makes the vertex and the momentum conservation factors 
complicate. 

(i) Draw a Feynman diagram. For an electron propagator, introduce the Green's function 



p U 

G { °1 = ► (C8) 



(ii) For one Coulomb line which is the interaction vertex, a local momentum interaction is not assigned in the usual 
way. We add a non-local interaction including four electron momentum with phase factor as 

PlUl \ k u 4Pi ^4 

V^k^e-^^-^y = )■■■{ ■ ( C9 ) 

Pi U 2 j Vp 3 ^3 

(iii) Add a momentum conservation factor for each interaction vertices with four electron momentum pi and 
frequency u>i 

^(27r) 2 (5( Pl - p 2 - k - 27rn)e^ (p2 ' n) (27r) 2 <5(p 3 - p 4 + k - 27rm)e l0(p3 ' m) 

nm 

x (2tt) 2 S(uj 1 ~uj 2 - uj)S(iv 3 - uj 4 + w). (CIO) 

The phase factor (p(pi,n) is added to the delta function. 

(iv) For one current vertex, a local momentum factor is also not assigned as the interaction vertex case. We add a 
non-local factor 

kuj 

g l ,{k)er^ k ' p ^ + ^ k ' k y = ► — X ( cu ) 

fl 

(v) Add a momentum conservation factor for each current vertices with two electron momentum p.; and frequency 

^(27r) 2 <S(pi - p 2 - k - 27rn)e l0(p2 ^ ) (27r)5(cj + wi - u 2 ). (C12) 

n 

The phase factor <fi(pi,n) is added to the delta function. 

(vi) Perform integral for internal momentum and add the numerical factor 

(ci3) 

where I is the number of internal line. Count the electron loop number and add the factor (— l) L for L electron loops. 
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APPENDIX D: DUALITY BETWEEN THE DIRECT TERM AND THE EXCHANGE TERM 

One of the surprising properties of the 2D electron system in a magnetic field is that the ladder diagrams take a 
similar form with the bubble diagram and the bubble and ladder diagrams can be summed up to the infinite order. 
This is caused by the duality between the direct term and the exchange term. In this section we show this property 
in the two-loop order as shown in Fig. ^] 

Following the rule of Appendix O the left two- loop diagram in Fig. ^| reads 



dcJidtJ2 / cPpi d 2 p 2 q(q) 



2vr 2vr J MBZ (2tt) 2 (2tt) 2 Pi-fe^i-w P2 -fc,^ 2 - w 



In general, the Fourier transform V;(q) = J (f^p ^Kp)e tpq satisfies the following duality relation between the direct 
term and the exchange term, 



Vi(p + 27rn)e i ( p+27rn '- q = ^(<1 + n)e- jpn . 



(D2) 



Using this relation, Eq. (|Dlf) is written as 
dwi duj^ r a2 ""- a2 



£ 



d 2 pi d 2 p 2 ^(0) 
2tt 2tt J mbz (2tt) 2 (2tt) 2 vx-k 

y^7r »(fc)W»(fc)7r ra0 (fc). 



2tt 



ria e 



>x(pi-p 2 ) 



(D3) 



This result is equivalent to two bubbles connected with the interaction W(k) as the right two-loop diagram in Fig. 1121 




Pi — k P2 - k 

-* = * <!>:-:<C> 4 

P\ P2 



FIG. 12: The two-loop ladder diagram and corresponding two-loop bubble diagram 
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